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SUMMARY 


A  method  is  presented  for  obtaining  amplitude  and  angle  plots  for 
rational  algebraic  functions  of  an  imaginary  variable.  Application  of 
the  method  is  illustrated  by  an  example  in  which  the  frequency  response 
of  an  automatically  controlled  aircraft  is  plotted.  The  method  involves 
the  use  of  templets  but  does  not  require  the  factoring  of  polynomials. 
Amplitude  and  angle  plots  of  high-degree  rational  functions  of  an 
imaginary  variable  can  be  obtained  more  rapidly  with  this  method  than 
by  analytical  calculation  or  by  methods  involving  the  factoring  of 
polynomials . 


INTRODUCTION 


If  a  dynamic  system  can  be  represented  by  a  set  of  linear  differ¬ 
ential  equations  with  constant  coefficients  and  with  time  as  the  inde¬ 
pendent  variable,  the  relationship  between  any  output  variable  of  the 
system  0o  and  any  input  of  a  system  9±  can  be  obtained  in  the  opera¬ 
tional  form 


=  K  bP  +  ap^eP'1  + 
sq  +  bg_^s q_1  + 


+  aO 

+  t>0 


where  s  denotes  the  operation  and  K,  aQ^  .  •  • ,  bQ^  •  •  •  > 

p,  and  q  are  constants.  This  function  is  usually  referred  to  as  the 
transfer  function  of  the  system.  It  is  shown  in  chapter  k  of  refer¬ 
ence  1  that  when  s  is  replaced  by  j<o,  where  j  =  fU  the  resulting 
complex  function  of  represents  the  frequency  response  of  the  system. 
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That  is ,  when  any  value  of  forcing  frequency  a),  in  radians  per  second, 
is  substituted  into  the  function,  the  amplitude  and  angle  of  the 
resulting  complex  value  of  the  function  represent  the  amplification 
and  phase  shift,  respectively,  through  the  system  when  excited  at  that 
frequency- 

In  the  analysis  and  synthesis  of  such  systems,  it  is  often  desira¬ 
ble  to  obtain  plots  of  the  amplification  or  amplitude  A  and  the  phase 
shift  or  phase  angle  <p  for  the  system  over  a  portion  of  the  frequency 
spectrum.  These  plots  can  be  obtained  by  evaluating  the  transfer  func¬ 
tion  of  jcs  at  discrete  values  of  cd  and  plotting  amplitude  and  angle 
as  functions  of  o>,  but  this  process  is  tedious  and  time  consuming. 
Another  method  of  obtaining  the  plots  is  to  factor  the  numerator  and 
denominator  of'  the  transfer  function  into  first-degree  and  second-degree 
factors  and  combine  them  graphically  by  using  templets  to  represent  each 
factor  on  logarithmic  coordinates  as  described  in  reference  1.  But, 
when  high-degree  polynomials  are  involved,  the  process  of  factoring 
becomes  very  tedious. 

The  present  method  allows  the  use  of  templets  for  plotting  rational 
functions  of  jo)  but  does  not  require  the  factoring  of  polynomials. 

The  basic  principle  involved  is  the  use  of  logarithms  (or  logarithmic 
graph  paper)  to  simplify  the  plotting  of  the  individual  terms  of  the 
real  and  imaginary  parts  of  the  numerator  and  denominator  polynomials. 
This  principle  was  initially  conceived  and  applied  to  rational  functions 
by  Mr.  William  M.  Kaufftaan  of  MCA  Headquarters.  Most  of  the  templets 
used,  in  the  present  method,  to  perform  the  graphical  operations  required 
in  obtaining  the  frequency  plots  were  first  developed  and  used  by 
Mr.  Kauffman. 


SYMBOLS 


ei 


3  -\£I 


output  variable  of  a  system 

input  (excitation)  of  a  system 

transform  variable  corresponding  to  differ¬ 
ential  operator  -sl¬ 
at 

excitation  frequency,  radians  per  second 


■0 
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= \j  |r(“3  d2  + 


=  Amplitude  of  denominator 


D 

^(“3  D 


=  Angle  of  denominator 


A(cd) 


N 

[a  (“3  D 


Amplitude  of  the  rational  function 


cp(co)  =  |5(“3n 


=  Angle  of  the  rational  function 


The  real  and  imaginary  parts  of  either  the  numerator  or  the  denomi¬ 
nator  can  be  obtained  by  algebraic  addition  of  the  individual  terms  each 

of  which  is  a  function  of  (D.  The  logarithm  of  a  typical  term  cafl  is 


log  co)n  =  log  c  +  n(log  cd) 


but  log  c  is  a  constant  and  n(log  co)  is  a  linear  function  of  log  co. 
Therefore ,  each  term  of  the  real  or  imaginary  part  of  numerator  or  denomi¬ 
nator  is  a  straight  line  of  slope  n  when  plotted  logarithmically  as 
a  function  of  cn,  that  is,  when  the  logarithm  of  the  function  is  plotted 
against  the  logarithm  of  co  on  a  linear  scale  or,  equivalently,  when  the 
function  is  plotted  against  co  on  a  logarithmic  scale.  (This  type  of 
plot  is  hereinafter  referred  to  as  a  logarithmic  plot. )  Although  a  loga¬ 
rithmic  plot  of  each  of  the  separate  terms  is  very  easy  to  make,  these 
terms  must  be  added  (or  subtracted),  and  graphical  addition  is  not  a 
simple  operation  when  the  curves  to  be  added  are  plotted  logarithmically. 
But  it  is  possible  to  make  a  templet  which  can  be  used  to  add  (or  subtract) 
graphically  any  two  values  plotted  logarithmically.  Such  a  templet  is 
shown  in  figure  1.  Its  construction  is  described  in  the  appendix  and  its 
practical  use  is  described  in  a  subsequent  section  entitled  "Practical 
Application  of  Methods  and  Description  of  Templets."  By  use  of  such  a 
templet,  any  two  logarithmic  curves  can  be  added  or  subtracted  point  by 
point. 
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Consider  any  two  functions  of  cjd,  F^(cjd)  and  F2(o>),  having 
values  P  and  Q,  respectively,  at  cd  =  and  plotted  logarithmically 
as  shown  in  figure  1.  Then 

ios[ji(“i) +  =  ios(p  +  q) 

■  log  q(i  ♦  £) 

=  log  Q  +  log  ^1  + 

Since  the  distance  between  P  and  Q  on  a  logarithmic  plot  depends 
only  on  the  ratio  P/q  and  not  on  the  individual  values  of  P  and  Q, 

a  templet  can  be  made  which  converts  this  distance  corresponding  to  log 

into  a  distance  corresponding  to  log  ^1  +  which  is  also  dependent  only 

on  the  ratio  P/Q,  and  adds  this  latter  distance  to  the  point  representing 
log  Q  to  locate  the  point  representing  log(P  +  Q).  Similar ly, 

log(P  -  Q)  =  log  P  +  log/l  -  ^  (P  >  Q) 


and  a  templet  can  be  made  which  converts  a  distance  corresponding  to 

log  |  to  a  distance  corresponding  to  log^l  -  and  adds  the  latter 

distance  to  the  point  representing  log  P  to  locate  the  point  repre¬ 
senting  log(P  -  Q).  Such  a  templet  is  shown  in  figure  1  and  is  herein¬ 
after  referred  to  as  the  summing  templet. 

If  the  logarithmic  plots  of  the  functions  F^co)  and  F2(cd)  are 

straight  lines,  as  in  the  case  of  the  individual  terms  of  the  real  and 
imaginary  parts  of  numerator  or  denominator  of  F(jco),  the  shape  of  the 
logarithmic  plot  of  their  sum  or  difference  can  he  shown  to  depend  only 
on  the  slopes  of  the  two  straight  lines.  That  is,  shifting  either  or 
both  of  the  straight  lines  in  any  direction  merely  shifts  their  point  of 
intersection  and  hence  shifts  the  curve  representing  their  sum  or  difference 
hut  does  not  alter  the  shape  of  the  curve.  Therefore,  templets  can  be 


P-tlOf* 
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made  representing  the  sum  and  difference  of  pairs  of  straight  lines  of 
various  slopes.  Each  of  these  templets  when  properly  located  can  give 
directly  a  logarithmic  plot  of  the  sum  or  difference  of  a  pair  of  terms 
of  the  real  or  imaginary  part  of  the  numerator  or  denominator  of  F(jo>). 
The  contours  for  four  such  templets,  which  are  sufficient  for  plotting 
polynomials  of  ninth  degree  or  less,  are  shown  in  figure  2  and  these 
templets  are  hereinafter  referred  to  as  term-pair  templets  or  term-pair 
contours . 

After  logarithmic  plots  of  the  real  and  imaginary  parts  of  numerator 
or  denominator  have  been  obtained,  it  is  necessary  to  obtain  the  corres¬ 
ponding  amplitude  and  angle.  The  amplitude  can  be  obtained  graphically 
by  use  of  a  templet  similar  to  the  summing  templet  since 


log  A(o>)  =  log  y[r(cd)]2  +  |i(cjd)J2 
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A  templet  can  be  made  which  at  any  value  of  frequency  converts  the 

I(“l)  /  R  (“]_)' 

'or  log  — ■ - — 4-1  into  a  distance 


distance  corresponding  to  log 


R 


K) 


K), 


corresponding  to  logi/1  + 


1 

1  + 

IK) 

2 

rK)_ 

or  log  1 / 1  + 


R(coi) 


IK) 


and  adds 


locate  the  point  representing  log  Such  a  templet 

is  shown  in  figure  3(a)  and  is  hereinafter  referred  to  as  the  amplitude 
templet. 


Since  the  angle  tan~^  is  dependent  only  upon  the  ratio 

R(cjo) 

I(cd)/r(o)),  it  can  he  obtained  graphically  at  any  value  of  frequency 
■j_  from  the  logarithmic  plots  of  I(cd)  and  R(cd)  by  use  of  a  scale 


CD 


which  converts  the  distance  corresponding  to  1 W  /RM  (  i.e, 9  the 
distance  between  i(cDp)  and  R(cdi)  on  a  logarithmic  plot)  into  angle 

directly.  Such  a  scale  is  shown  in  figure  3(b)  and  is  hereinafter  referred 
to  as  the  angle  scale 

After  obtaining  the  amplitudes  Ja(cdJJ^  and  |a(cd7|  ^  and  the 
angles  £q)((D)J  N  311(1  (9(ai)J  j),  it  is  then  necessary  to  combine  them  to 

obtain  the  amplitude  A(cd)  and  the  angle  cp(cd)  of  the  original  function 

F(joo)  =  KA(co)e^  m  „  Nc  ^  J 


The  angle  is  obtained  by  merely  subtracting  |9(cd)J  D  ^rom  [<p(c)j]  ^ 
for  any  value  of  frequency  cDj_  after  [t(cd7|  p  and  N  have  "been 

obtained  by  use  of  the  angle  scale. 


By  taking  logarithms  of  both  sides  of  the  relation  A(cd)  = 
the  following  relation  is  determined: 


N 


[A(mj]  D 


log  A  (to)  =  log  jA(<x.Jj  n  "  log[A(a>yjD  =  logjA(m)j  N  +  log— ■■■— 
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Therefore ,  a  logarithmic  plot  of  A(<x>)  can  he  obtained  by  graphically 
subtracting,  point- by- point,  the  logarithmic  plots  of  jj  and  P<  0) 


or  "by  adding  the  logarithmic  plots  of 


Since  log  1=0,  the  axis  representing  unity  on  the  logarithmic 
plots  must  be  used  as  the  reference  or  zero  axis  when  adding  or  subtracting 
distances  on  a  logarithmic  plot.  The  logarithmic  plot  of  l/f(o>)  is 
merely  the  logarithmic  plot  of  f(o>)  reflected  in  the  axis  representing 
unity. 


PRACTICAL  APPLICATION  OF  METHOD  AND  DESCRIPTION  OF  TEMPLETS 


The  general  form  of  either  the  numerator  or  denominator  polynomial 


of  a  rational  algebraic  function  of  (jo)),  F(Jcd)  = 


P2W 


is 


p(j^) 


a0  +  al('«H  +  a2(jo))2  +  a3(jao)3  +  a^(joa)^  + 

a^Cj cd)5  +  ag(jo))6  +  a„(>)7  +  .  .  . 

ao  +  Oa]5>  -  a2<n2  -  ja^o^  +  +  ja^o)5  _  agm^ 


where  aQ,a-^,a2,  .  .  are .  constants .  In  the  practical  case,  the 

variable  03  is  usually  excitation  frequency.  In  certain  cases,  it  may 
be  necessary  to  factor  (joo)11  from  the  polynomial  to  obtain  a  constant 
term  in  the  remaining  factor  as  in  the  above  form  and  to  correct  later  for 
the  factor  (jcs)n.  This  point  is  mentioned  again  subsequently. 

The  real  and  imaginary  parts  of  P(jco)  are 

Real  part  =  R(cjd)  =  ^aQ  -  a2a^')  +  (^k0^  “  a6a3^)  +  •  •  • 

Imaginary  part  =  l(o o)  =  -  a^oy^  +  +  .  .  . 
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and  term-pair  templets  can  be  made  to  represent  each  of  the  term-pairs 
enclosed  by  parentheses.  The  contours  of  such  templets  are  shown  in 
figure  2. 

In  order  to  get  either  R(tu)  or  I (00) ^  the  proper  term-pairs  are 
plotted  by  use  of  the  term-pair  templets  and  are  summed  by  use  of  the 
summing  templet.  In  order  to  use  the  term-pair  templets,  the  inter¬ 
section  of  the  straight  lines  representing  the  logarithmic  plots  of  the 
separate  terms  ±s  first  located.  This  intersection  fixes  the  position 
of  the  term-pair  templet  so  that  its  contours  represent  the  logarithmic 
plot  of  the  sum  and  difference  of  the  pair  of  terms.  If  ao,ap,a2,  .  .  . 

are  all  of  the  same  algebraic  sign,  the  terms  of  R(u>)  and  I  (to)  will 
have  alternating  algebraic  signs  so  that  term-pairs  composed  of  consecu¬ 
tive  terms  will  always  be  the  difference  of  two  terms  rather  than  the 
sum. 


It  may  be  noted  here  that  logarithmic  plots  can  represent  only 
absolute  (or  positive)  values.  It  is  apparent  that  for  large  values  of 
co,  the  expression  aQ  -  a2oy^  for  example,  will  be  negative  (assuming  aQ 

and  a2  to  be  positive).  For  these  values  of  co,  a  logarithmic  plot 

must  represent  agco^  -  aQ  or  -  |a.Q  -  instead  of  aQ  - 

When  graphically  obtaining  the  sum  or  difference  of  two  functions  F1(co) 
and  F2(co),  it  is  necessary  to  know  whether  the  plot  represents  F(co) 

or  -F (co) .  If  then,  for  example,  a  logarithmic  plot  of  Fp(co)  +  F2(co) 

is  to  be  obtained  and  logarithmic  plots  of  F1(co)  and  -F2(co)  have  been 
obtained,  Fp(co)  +  F2(co)  may  be  found  by  graphically  subtracting  -F2(co) 

from  F-]_(co)  since  F^co)  +  F2(co)  =  F^(co)  -  £-F2(co)J  . 

After  two  logarithmic  curves  (such  as  term-pair  contours)  are  drawn, 
their  sum  or  difference  can  be  obtained  graphically  by  use  of  the  summing 
templet  of  figure  1.  This  templet  is  used  to  add  two  logarithmic  curves 
F-^(co)  and  F2(gd)  graphically  as  follows: 

The  reference  point  X  (a  small  hole  in  the  templet)  is  set  on  the 
lower  curve  at  some  value  of  frequency  o)p,  and  contour  a  is  set  on 

the  upper  curve  at  CD]_  as  illustrated  in  figure  1.  The  point  representing 

the  sum  of  the  two  curves  at  a)]_  is  found  on  contour  b  at  a>]_.  This 

process  is  repeated  at  as  many  values  of  co  as  are  needed  to  define  the 
resulting  logarithmic  plot. 
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In  order  to  subtract  two  logarithmic  curves,  the  reference  point  X 
is  set  on  the  upper  curve  at  some  value  of  frequency  goj_  and  contour  c 

(or  contour  d)  is  set  on  the  lower  curve  at  00-^.  The  point  representing 
the  difference  of  the  two  curves  at  GOp  is  found  on  contour  d  (or 
contour  c)  at  GOp.  This  point  represents  the  upper  curve  minus  the 

lower  curve  or  the  negative  of  the  lower  curve  minus  the  upper  curve 
at  GDp.  Repeat  at  as  many  values  of  go  as  are  needed  to  define  the 
resulting  logarithmic  plot. 

This  process  of  addition  and  subtraction  of  curves  is  continued 
until  all  terms  of  R(go)  or  I  (go)  have  been  included  and  logarithmic 
plots  of  R(go)  and  I  (go)  have  thus  been  obtained.  Then 

+  [r(go)|  2  can  be  obtained  graphically  by  applying  the  amplitude 

templet  of  figure  3(a)  to  the  logarithmic  plots  of  R(cjo)  and  I (cjo)  as 
follows: 

The  reference  point  Y  of  the  amplitude  templet  (represented  by  a 
small  hole  in  the  templet)  is  set  on  the  lower  curve  at  some  frequency  GOp, 

and  contour  e  is  set  on  the  upper  curve  at  GOp.  The  point  representing 

at  go-^  is  on  contour  f  at  go^.  This  process  is 
repeated  at  as  many  values  of  go  as  are  needed  to  define  the  resulting 

logarithmic  plot  of 


The  angle  scale  of  figure  3(b)  is  used  in  conjunction  with  logarithmic 

plots  of  R(go)  and  I  (go)  to  obtain  values  of  tan”1  ^4^4*.  The  line  r 

R(go) 

near  the  middle  of  the  scale  is  placed  on  the  logarithmic  plot  of  R(oo) 

-1  I  fa) 

at  any  frequency  value  go-i  .  The  angle  tan  x  is  read  directly 

Eh)  .  . 

in  degrees  from  the  scale  at  the  point  representing  I(GOp).  The  quadrant 

of  the  angle  is  determined  from  the  algebraic  signs  of  I  and  R  as 
indicated  on  the  angle  scale.  This  process  is  repeated  for  as  many  values 
of  go  as  are  needed.  The  angle  values'  thus  obtained  may  be  either  tabu¬ 
lated  or  plotted  as  a  function  of  go. 


is  read  directly 
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The  remaining  step  is  the  combination  of  numerator  and.  denominator 
to  obtain  the  function 


F(Jcd)  =  KA(co)ejcP(“) 

[aw]  j{Ij(“3h  -  M,} 

°kmt 

The  angle  cp(m)  may  be  obtained  by  point-by-point  subtraction 
of  [cpM]])  from  [cp(“)Jw  either  graphically  or  analytically.  The 
logarithmic  plot  of  A(u>)  may  be  obtained  by  point-by-point  graphical 
addition  of  log[A(4]N  and  log  — (or  logfl (coTJ N  -  log[A(a>y| D) . 

If^  (*>)*  was  factored  from  either  the  numerator  or  denominator 
polynomial,  it  must  be  included  at  this  point.  That  is,  the  amplitude 

and  the  angle,  n  X  90°  must  be  added  to  or  subtracted  from  the 
amplitude  and  angle  curves  already  obtained. 

For  a  polynomial  P(ja>)  =  aQ  +  a^jco)  +  a2(jo>)2  +  .  .  the  values 

of  the  coefficients  ao>al->a2  .  .  .  may  be  of  such  widely  differing 

magnitudes  as  to  make  the  plotting  of  the  various  terms  on  the  same  graph 
impractical.  This  situation  may  be  remedied  by  a  change  in  frequency 
scale  to  =  ku.  This  and  the  other  foregoing  procedures  are  illustrated 
m  the  section  entitled  "Illustrative  Example." 


Illustrative  Example 


A  transfer  function  representative  of  those 
systems  is 


encountered  in  physical 


F(s)  = 


_ 73l(s3  +  91.9s2  +  8o60s  +  26850) _ 

+  93. 5s^  +  9HOs3  +  273000s2  +  7550000s  +  19630000 
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This  transfer  function  was  encountered  in  analyzing  an  idealized  system 
made  up  of  an  airframe  and  automatic  control  system  which  controlled  the 
pitch  attitude  of  the  airframe. 


If  s  is  replaced  by  joo  where  to  is  excitation  frequency  and 

J  is  \Ti,  the  resulting  function  of  jo)  represents  the  steady— state 
frequency  response  of  the  system  (reference  1,  p.  96) 


F(jcd)  = 


731  |Tjcu) 3  +  91- 5  ( +  8060  (jcn)  +  26850] 


(jo>)5  +  93.5(jtD)^  +  9110  (jo>)  3  +  273000  ( joo) 2  +  7550000  (jco)  +  19630000 


The  magnitudes  of  the  coefficients  of  the  separate  terms  differ  greatly, 
the  ratio  of  the  largest  and  smallest  being  19,630,000.  This  makes  the 
plotting  of  the  various  terms  on  the  same  graph  inconvenient.  However, 
if  a  change  in  frequency  scale  o>  =  40u  is  made,  the  resulting  expression 


,  ,  0.457  (ju)3  +  2.29(ju)2  +  5.04(Ju)  +  O.42] 

F(ju)  =  - 4= - _J _ 

(ju)5  +  2.34(ju)1<-  +  5-  69(  ju)3  +  4.27(ju)2  +  2.95(ju)  +  0.192 
The  ratio  of  the  largest  and  smallest  coefficients  in  this  expression 

5.69  . 

1S  0  192  =  The  factor  °*lt-57  being  ignored  for  the  present,  the 

real  and  imaginary  parts  of  numerator  and  denominator  are 


[?(Ul]  N  =  (°-^2  “  2-29u2) 

|j(uT]N  =  (5. 04u  -  u3) 

|r(uI]d  =  (0.192  -  4.27u2)  +  2. 34u^ 

[i(uy]D  =  (2  .  95u  -  5.69u3)  +  u5 
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The  contours  of  the  two  templets  of  figures  2(a)  and  2(h)  can  he 
used  to  represent  the  term-pairs  in  parentheses.  Since  the  numerator 
is  only  of  third  degree,  its  real  and  imaginary  parts  can  each  be  repre¬ 
sented  directly  by  one  of  the  term-pair  templet  contours.  These  contours 
are  shown  properly  Orientated  on  the  graph  in  figure  4,  and  the  real 
and  imaginary  parts  are  combined  by  use  of  the  amplitude  templet  to  obtain 

[A(u)]N  =  \J [i (u)j  jj2  +  |R(ugN2  at  any  value  of  u.  Figure  4  illustrates 


is  obtained  from  the 


the  use  of  the  amplitude  templet  to  obtain  _V(0.365)2  +  (0.755)2  =  0.84 
at  u  =  0.15.  The  angle  f9(ujlw  =  tan""*-  — — =1^  is  obtained  from  the 

[r(u3  n 

logarithmic  plots  of  the  real  and  imaginary  parts  by  use  of  the  angle 
scale.  Figure  4  illustrates  the  use  of  the  angle  scale  to  obtain 

tan  ^  ~ ~~ v?  =  114.5°  at  u  =  1.  The  angle  values  are  tabulated  for  dis- 
-1. 84 

Crete  values  of  u  in  table  I. 

Since  the  denominator  of  F(ju)  is  of  fifth  degree,  an  extra  term 
must  be  added  to  a  term-pair  to  obtain  either  the  real  or  imaginary  part 

of  the  denominator.  This  addition  is  performed  in  figure  5  for  ^(u^j^ 

and  in  figure  6  for  £l(uy|  D*  Figure  5  illustrates  the  use  of  the  summing 

templet  for  graphically  subtracting  2.5  -  0. 9  =  1.6.  Figure  6  illustrates 
the  use  of  the  summing  templet  for  graphically  adding  0.120  +  0.335  =  0.455 
The  amplitude  templet  and  the  angle  scale  are  then  applied  to  the  loga¬ 
rithmic  plots  of  [R(u)]d  and  Jl(u)]  D  to  obtain  [a(uTJ  d  and  (rUTjj) 

shown  in  figure  7  and  table  I,  respectively.  The  operations  of  fig¬ 
ures  4,  5 ,  6,  and  7  could  be  carried  out  on  a  single  graph,  but  they  have 

been  separated  here  for  clarity.  Values  of  £cp(u)|  ^  are  tabulated  in 
table  I  for  discrete  values  of  u. 

The  subtraction  of  Pu3d  from  |r(uT]  n  yields  total  angle  cp(u) 
which  is  also  tabulated  in  table  I.  If  |A(ujlN  and  - i -  are 

L  J  (A(u7]D 

plotted  on  the  same  graph,  they  can  be  combined  to  obtain  A(u)  as  in 
figure  8.  This  operation  consists  of  the  addition  of  the  graphical 
distances  measured  vertically  from  the  horizontal  axis  representing  unity 


which  is  also  tabulated  in  table  I.  If 
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(since 


log  1  =  0)  to  the  logarithmic  plots  of 


Distances  measured  upward  from  the  unity  axis  are  considered  positive 
while  those  measured  downward  from  the  unity  axis  are  considered  negative. 
The  factor  0.457  of  F(ju),  which  has  been  ignored  until  now,  has  been 
taken  into  account  at  this  point  by  shifting  the  logarithmic  plot  of 


vertically  an  amount  corresponding  to 


log  0.457. 


Plots  of 


A(u)  and  cp(u)  obtained  by  the  present  method  are  compared  with  the 
analytically  calculated  plots  in  figures  8  and  9* 


General  Comments 


Because  of  the  steep  slopes  of  the  curves  being  added  and  subtracted 
and  because  this  method  frequently  involves  the  graphical  subtraction  of 
two  large,  nearly  equal  quantities,  care  in  performing  the  graphical 
operations  is  necessary  in  order  to  maintain  a  high  degree  of  accuracy. 
After  some  familiarity  with  this  method  is  acquired,  however,  amplitude 
and  angle  plots  for  high-degree  rational  functions  of  jo)  can  be  obtained 
more  rapidly  with  this  method  than  by  analytical  calculation  or  by  fac¬ 
toring  the  polynomials  and  plotting  with  first-  and  second-degree  templets. 
The  inaccuracies  usually  occur  in  the  upper  frequency  range.  If  such 
inaccuracies  can  be  tolerated,  the  time  required  to  obtain  plots  can  be 
decreased  considerably. 


One  possible  method  of  decreasing  the  slopes  of  the  individual  terms 
of  the  real  and  imaginary  parts  of  an  nth  degree  polynomial  might  be  to 


factor 


( joo)m  from  the  polynomial  where 


This  process  results  in 


a  modified  polynomial  whose  real  and  imaginary  parts  contain  terms  of  both 
positive  and  negative  slopes  when  plotted  logarithmically.  After  the 
amplitude  and  angle  of  this  modified  polynomial  have  been  obtained,  they 

must  be  corrected  by  the  factor  (ju>)m  to  obtain  the  amplitude  and  angle 
of  the  original  polynomial.  Obtaining  amplitude  and  phase  of  the  modified 
polynomial  by  the  method  mentioned  here  would  require  term-pair  templets 
representing  terms  having  both  positive  and  negative  slopes.  The  possible 
advantages  of  this  refinement  have  not  been  thoroughly  investigated. 


The  use  of  logarithmic  graph  paper  different  from  that  used  here  may 
be  desirable  for  many  problems.  For  most  problems,  the  use  of  logarithmic 
paper  having  more  cycles  in  the  vertical  direction  has  been  found  advanta¬ 
geous.  Templets  are  easily  constructed  for  use  with  any  logarithmic  scale 
chosen  (see  appendix). 


* 


1 6 


MCA  TN  2522 


CONCLUDING  REMARKS 


A  method  has  been  presented  for  obtaining  amplitude  and  phase  plots 
for  rational  algebraic  functions  of  an  imaginary  variable  F(joi).  The 
method  involves  the  use  of  templets  but  does  not  require  the  factoring 
of  polynomials.  After  some  familiarity  with  the  method  is  acquired, 
amplitude  and  angle  plots  of  high-degree  rational  functions  of  jcn  can 
be  obtained  more  ^apidly  with  this  method  than  by  analytical  calculation 
or  by  methods  that  involve  the  factoring  of  polynomials. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  July  24,  1951 


* 


MCA  TN  2522 


17 


APPENDIX 

FASHIONING  THE  TEMPLETS 


The  summing  templet,  the  amplitude  templet,  and  the  angle  scale  each 
makes  graphical  use  of  two  quantities  q^_  and  qg  to  obtain  a  desired 

result.  It  has  been  shown  in  the  section  entitled  "Theoretical  Analysis" 
that  the  use  of  either  of  these  templets  is  independent  of  the  magnitude 
of  the  individual  quantities  concerned  but  depends  only  on  their  ratio, 
that  is,  for  example,  the  same  points  on  the  summing  templet  can  be  used 


to  obtain  1+2=3  or  3+6=9 


^s ince 


Therefore,  if  one 


of  the  quantities  q^  is  arbitrarily  chosen  to  be  unity,* the  values  of 
the  result  (q^  +  q2  in  the  case  of  the  summing  templet)  can  be  calculated 
and  tabulated  for  discrete  values  of  the  second  quantity  q2.  The  values 

so  tabulated  can  then  be  used  with  the  chosen  logarithmic  scale  to  make 
the  templet.  A  tabulation  of  the  values  of  1  ±  q2  for  use  in  making 

the  summing  templet  is  not  necessary  since  these  values  are  easily  calcu¬ 


lated  mentally.  A  tabulation  of  the  values  of  f  +  (q2)2  fo  r  use  in 
making  the  amplitude  templet  is  given  in  table  II,  and  a  tabulation  of 


the  values  of  q2  against  tan' 
is  given  in  table  III. 


1  q2 


for  use  in  making  the  angle  scale 


The  exact  shape  of  the  contours  of  the  summing  templet  or  the  ampli¬ 
tude  templet  is  arbitrary,  the  only  requirement  being  that  the  desired 


relationship  between  q^,  q^ 
satisfied  for  all  values  of 


q2 


and  qjL  ±  qg 


(or  \JN2  +  N2) 


be 


by  the  contours  of  the  templets  when 


applied  to  the  chosen  logarithmic  scale.  The  general  shape  chosen  here 
is  convenient  to  use  because  the  reference  point  can  be  anchored  to  a 
desired  point  on  the  graph  with  a  sharp  pointed  instrument,  and  the 
templet  can  be  rotated  about  this  point  to  the  desired  position. 


The  term-pair  templets  all  represent  frequency  functions  of  the 

form  ±aaP‘  -  bo>n+2>  a s  was  pointed  out  in  the  section  entitled  "Theo¬ 
retical  Analysis,"  the  shape  of  the  logarithmic  plot  of  such  a  function 
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is  independent  of  the  magnitudes  of  a  and  "b.  Therefore,  the  shape 
of  the  contour  representing  such  a  function  plotted  to  the  chosen 
logarithmic  scale  can  he  determined  by  assuming  a  =  b  =  1.  Values  of 

-an)n  ^  bo)n+^  can  then  be  tabulated  for  discrete  values  of  for  any 
desired  value  of  n.  Such  a  tabulation  is  given  in  table  IV  for 
n  =  0,  1,  k,  and  5.  The  values  tabulated  are  sufficient  to  make  the 
four  term-pair  templets  required  for  plotting  polynomials  of  ninth  degree. 
By  plotting  these  values  to  the  chosen  logarithmic  scale ,  the  contours 
of  the  term-pair  templets  are  obtained.  These  values  are  plotted  in 
figure  2  to  the  logarithmic  scale  used  in  this  report  which  is  the  scale 
of  Keuffel  &  Esser  Co.  No.  359-HlCr  logarithmic  graph  paper. 

Satisfactory  templets  can  be  made  from  thin,  transparent  sheet 
plastic  which  can  be  cut  with  scissors  and  smoothed  with  a  fine  file. 
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TABLE  I 


PHASE  ANGLE  VALUES  FOR  ILLUSTRATIVE  EXAMPLE 


u 

(40  radians/sec) 

Hn 

(deg) 

(deg) 

cp(u) 

(deg) 

<p(u) 

(deg) 

(a) 

0.07 

42 

51 

-9 

-9.2 

.08 

46 

55 

-9 

-9.7 

.09 

49 

58.5 

-9-5 

-10.2 

.1 

52 

62.5 

-10.5 

-10.8 

•  15 

64 

77 

-13 

-13.8 

.2 

71.5 

87 

-15.5 

-15.5 

.3 

82 

103 

-21 

-21.5 

.4 

— 

-29.3 

•  5 

93 

132 

-39 

-39-2 

.  6 

98 

147.5 

-49.5 

-51.6 

.7 

102.5 

166.5 

-64 

-66.1 

.8 

106.5 

186 

-80.5 

-81.6 

•  9 

111 

206 

-95 

-96.3 

1.0 

114.5 

225 

-110.5 

-108. 8 

1.1 

119 

239 

-120 

-119.0 

1.2 

123 

252 

-129 

-127. 3 

1.3 

128 

265 

-137 

-134.3 

1.4 

132.5 

277-5 

-145 

-140.5 

1.5 

138 

289 

-151 

-146.2 

Calculated  analytically. 
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TABLE  II 

VALUES  FOR  USE  IN  MAKING  THE  AMPLITUDE  TEMPLET 


q2 

q2 

V1  +  (<l2) 2 

sj1  +  (^2) 2 

1.0 

1.414 

5.6 

5.689 

1.2 

1.562 

5-8 

5.886 

1.4 

1.T20 

6.0 

6.083 

1.6 

1.887 

6.2 

6.280 

1.8 

2.059 

6.4 

6.478 

2.0 

2.236 

6.6 

6.675 

2.2 

2.417 

6.8 

6.873 

2.4 

2.600 

7.0 

7.071 

2.6 

2.786 

7-2 

7.269 

2.8 

2.973 

7.4 

7.467 

3.0 

3.162 

7.6 

7.655 

3-2 

3-.  353 

7.8 

7.864 

3.4 

3.544 

8.0 

8.062 

3.6 

3.736 

8.2 

8.261 

3.8 

3.829 

8.4 

8.459 

4.0 

4.123 

8.6 

8.658 

4.2 

4.317 

8.8 

8.857 

4.4 

4.512 

9.0 

9.055 

4.6 

4. 707 

9.2 

9.254 

4.8 

4.903 

9.4 

9.453 

5.0 

5.099 

9.6 

9.652 

5.2 

5.295 

9.8 

9.851 

5.4 

5.492 

10.0 

10.05 
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TABLE  III 

VALUES  FOR  USE  IN  MAKING  THE  ANGLE  SCALE 


tan'^qg 

(deg) 

q2 

tan- ^2 
(deg) 

q2 

tan“-*-q2 

(deg) 

12 

1 

0.0175 

31 

0.6009 

61 

1.804 

2 

.0350 

32 

.6249 

62 

1.881 

3 

.0524 

33 

.6494 

63 

1.963 

4 

.0699 

34 

.6745 

64 

2.050 

5 

.0875 

35 

.7002 

65 

2.144 

6 

.1051 

36 

.7265 

66 

2.246 

7 

.1228 

37 

.7536 

67 

2.356 

8 

.1405 

38 

•  7813 

68 

2.475 

9 

.1584 

39 

.8098 

69 

2.605 

10 

.1763 

4o 

.8391 

70 

2.748 

11 

.1944 

4l 

.8693 

71 

2.904 

12 

.2126 

42 

.9004 

72 

3.078 

13 

.2309 

43 

.9325 

73 

3.271 

14 

.2493 

44 

.9657 

74 

3.487 

15 

.2680 

45 

1.000 

75 

3.732 

16 

.2868 

46 

1.036 

76 

4.011 

17 

•  3057 

47 

1.072 

77 

4.332 

18 

•  3249 

48 

1.111 

78 

4.705 

19 

.3443 

>9 

1.150 

79 

5.145 

20 

.3640 

50 

1.192 

80 

5.671 

21 

.3839 

51 

1.235 

81 

6.314 

22 

.4040 

52 

1.280 

82 

7.115 

23 

.4245 

53 

1.327 

83 

8.144 

24 

.4452 

54 

1.376 

84 

9-514 

25 

.  4663 

55 

1.428 

85 

11.43 

26 

.4877 

56 

1.483 

86 

14.30 

27 

.5095 

57 

1.540 

87 

19.08 

28 

.5317 

58 

1.600 

88 

28.64 

29 

•  5543 

59 

1.664 

89 

57.29 

30 

•  5774 

60 

1.732 

90 

0 
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VALUES  FOR  USE  IN  MAKING  TERM-PAIR  TEMPLETS 
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Figure  1.-  Summing  templet  with  illustration  of  its  use 
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(b)  Contours  of  term-pair  templets  representing  logarithmic  plots  of 
±au  ±  bu3^  shown  oriented  for  a  =  b  =  1. 

Figure  2*-  Continued* 
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(c)  Contours  of  term-pair  templets  representing  logarithmic  plots  of 
±au^  ±  bu^,  shown  oriented  for  a  =  b  =  1. 

Figure  2.-  Continued. 
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(d)  Contours  of  term-pair  templets  representing  logarithmic  plots  of 
±au^  -  hvJ f  shown  oriented  for  a  =  b  =  1. 

Figure  2.-  Concluded. 


Figure  3.-  Amplitude  templet  and  angle  scale 
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u,  k-radians  per  second 


Figure  U.-  Logarithmic  plots  of  -[£(u)]lP  ±[J(uI]n*  ^or 

illustrative  example  showing  the  use  of  the  amplitude  templet  and  the 
angle  scale. 
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Figure  Logarithmic  plots  of  ±(0.192  -  ii.27u2),  2.3W+,  and 

for  illustrative  example  showing  the  use  of  the  summing  templet  for 
subtraction. 
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Figure  6.-  Logarithmic  plots  of  ±(2.95>u  -  £.69u3),  u^,  and  -[i(u3d 

for  illustrative  example  showing  the  use  of  the  summing  templet  for 
addition. 
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Figure  7.-  Logarithmic  plots  of  ± |r!( uJJ j  [a(u2d 

illustrative  example* 
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Figure  Q.-  Plot  of  cp(u)  for  illustrative  example. 
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